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Abstract 

We prove the asymptotic stability of 1:1 spin orbit resonance for a 
stiff viscoelastic rotationally invariant satellite. The main point is that 
we do not assume any dissipation on the orbital degrees of freedom and 
we prove that anyway they relax to those of a circular orbit, due to the 
friction acting on the internal degrees of freedom of the satellite. Tech- 
nically the result is obtained by using the principal moments of inertia 
as coordinates in the space of elastic configurations and in proving the 
asymptotic stability through LaSalle's principle, using the energy as a 
Lyapunov function. 



1 Introduction 

From the point of view of dynamical systems, the dynamics of a viscoelastic 
satellite interacting with a pointlike planet should be studied by writing down 
the equation of motion for the internal degrees of freedom of the satellite (which 
in general is a partial differential equation) , completing the system by adding the 
equations for the orbital degrees of freedom and describing the corresponding 
solutions. Spin orbit resonance should appear as an asymptotic state of the 
system when the internal degrees of freedom are subject to some dissipation. 
We emphasize that the orbital degrees of freedom must not be subject to any 
dissipative force, since the orbital motion of the satellite is not subject to any 
direct friction. Relaxation of the orbital degrees of freedom should be proved to 
appear as a consequence of the interaction between internal and orbital degrees 
of freedom. 

In the present paper, we consider a general model of viscoelastic satellite in- 
teracting in quadrupole approximation with the gravitational field of the planet. 
We assume that the satellite is rotationally invariant, possesses a nondegenerate 
equilibrium configuration (under elastic and selfgravitating forces) and we study 
dynamics in the limit of large stiffness. Furthermore we restrict to the "planar 
situation" in which the center of mass lies in a plane, the spin axis is orthogonal 
to such a plane and the deformations are such that one of the principal axes of 
inertia of the body is always orthogonal to the plane of the orbit. We assume 
that the elastic degrees of freedom are subject to a viscoelastic force and we 
prove that the circular orbit in which the spin exhibits 1:1 resonance with the 
orbital frequency is asymptotically stable in the sense of dynamical systems. No- 
tice that also the orbital degrees of freedom are subject to relaxation, so that, in 
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particular, the eccentricity of the orbit tends to zero. In plane words, we prove 
that, if the initial data are sufficiently close to those corresponding to a circular 
orbit with 1:1 resonance, then the system relaxes to such an asymptotic state. 
However, our result is not in contradiction with the fact that, with arbitrary 
initial data, the dynamics may exhibit a different behaviour; furthermore, we 
do not have any estimate on the time needed to observe relaxation. 

The starting point of the proof is the remark that, in quadrupole approx- 
imation, the gravitational potential of an extended body in an external gravi- 
tational field is a function only of the principal moments of inertia and of the 
directions of the principal axes of inertia. So it is natural to use such quantities 
as coordinates in the configuration space of the elastic satellite. Furthermore, 
a good coordinate system has also to distinguish between the internal degrees 
of freedom (which are subject to friction) and the external ones (rotations and 
translations), which are NOT subject to friction. Technically the introduction of 
a good coordinate system is the most difficult part of the proof. The difficulties 
one has to face are of three kinds: the first one, which will be discussed in detail 
in Section [51 is that the separation between the rotational degrees of freedom 
and the elastic ones relies on an arbitrary choice, so one has to show that the 
final result does not depend on the particular choice made. The second difficulty 
is related to the fact that the principal axes of inertia do not define directly a 
coordinate system, since they are defined up to orientation. A third difficulty is 
that the system of coordinates we are aiming to introduce is singular at those 
configurations where the three principal moments of inertia are not distinct. 

As a first result we show that the principal moments of inertia and a couple 
of matrix of rotation related to the direction of the principal axes of inertia can 
be completed to a system of coordinates on the configuration space of the body, 
after removing the set of singular configurations (see Sect. 2). 

Then we have to write down the Lagrangian of the system in these coordi- 
nates. To this end we exploit the rotational invariance of the body: surprisingly 
enough this allows one to show that the elastic potential has a very special 
form which can be written down explicitly, and also to deduce some important 
properties of the kinetic energy. 

The subsequent steps are quite standard. First we consider the conservative 
dynamics of the system: we write the Lagrangian, exploit the conservation of 
angular momentum in order to reduce by one the number of degrees of freedom 
and prove that the family of solutions corresponding to 1:1 resonance forms a 
Lyapunov stable manifold of the reduced system. Finally we add friction on 
the internal degrees of freedom and prove that such a manifold becomes asymp- 
totically stable. The proof of the asymptotical stability is done by means of 
LaSalle's invariance principle, which is a refinement of Lyapunov's theory. The 
method substantially consists in observing that the energy is a nonincreasing 
function of time and showing that no orbit, except for the points of the mani- 
fold corresponding to 1:1 resonance, is contained in the submanifold of the phase 
space where dissipation does not occur. 

It is worth mentioning that, even if we prove that the speed of rotation of the 
satellite tends to zero, we are unable to prove that asymptotically "the satellite 
stops rotating": we only prove that asymptotically the principal axes of inertia 
stop rotating. This is compatible with a situation in which the satellite has a 
constant shape, but the direction of the deformation continuously changes in 
the body. In a pictorial way one can think of a rubber baloon which slides more 



2 



and more slowly on a wooden cross whose axes are fixed in space. 

We conclude this introduction discussing the connection with some previous 
results on spin orbit resonance, starting from the pioneering work by Darwin 
(see also [7]). Darwin's theory is based on the idea of phase lag, a quantity 
which describes the delay between the phase of tidal waves and the phase related 
to the change of direction of the axes of deformation of the body. Such a phase 
lag incorporates all the effects of friction. The methods and the ideas used in 
order to estimate the phase lag lead to different theories which were developed 
by many authors [TllSl fTUlfT^fH] (see also [7]). Sometimes the ideas of Darwin's 
theory are used to introduce an effective dissipative force on the spin degrees 
of freedom (see also [H]). Our approach is much more fundamental, since 
we actually study the equations of motion of the viscoelastic satellite, without 
introducing any intermediate concept as the phase lag. Since our theory is 
based on first principles, it could be interesting to develop it further in order to 
possibly justify Darwin's theory, at least in some particular situations. 

A series of studies taking a point of view closer to ours have been devel- 
oped in the papers [T21[TB]. In particular, in [T5] the dynamics of a homo- 
geneous viscoelastic satellite in an external gravitational field is studied in a 
setting very similar to ours. In such papers results are also deduced within a 
dynamical system perspective, and there is even no restriction to the "planar 
situation". Nevertheless, several approximations are made in the study of the 
evolution equations, in particular the terms involving second time derivatives 
are neglected, an operation whose possible rigorous justification is very delicate. 
Furthermore the analysis is performed using the averaging principle in a formal 
way. 

We also recall the papers , where the spin orbit problem is studied from 
the point of view of dynamical systems, getting remarkable results. However, 
they do not consider the same problem as in the present paper: they always put 
the satellite on a fixed elliptic orbit, and moreover they consider it as a rigid 
body, possibly subject to some effective dissipative forces. 

Acknowledgements. We thank Antonio Giorgilli, Alessandra Celletti and Sylvio 
Ferraz-Mello for some discussion on this problem, which led to considerable 
improvements of the result. 

2 Coordinates on the configuration space of an 
elastic body 

In this section we study the kinematics of an elastic body. The aim is to intro- 
duce suitable coordinates for the study of our problem. 

2.1 General considerations 

We will use the Lagrangian description of elasticity. In this approach one defines 
the so called material space 57, which is essentially an abstract realization of the 
elastic body in some reference configuration. In our case we choose il to be 
a three dimensional sphere. In the following we will assume that the body 
is invariant under rotations, in particular we assume that its density function 
p : il ^ is invariant under rotations. We denote by m the mass of the elastic 
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body, i.e. 

m = p{x)(fix . 
Jn 

The configuration of tlie body is a map C : J7 ^ , wliich gives the position 
in space of the the point x £ fi. 

Of course the map ( describes both the deformation of the body and its 
position and rotation in space, so it is natural to try to decompose ( into a rigid 
translation, a rigid rotation and an internal deformation. To start, define the 
center of mass of the body by 

X = 1 / C(x)p(x)d3x , (2.1) 
m Jn 

and decompose the configuration vector field <^ as 

C(x) - X + z;(x) , (2.2) 

where v is such that 

/ v{x)p{yi)(fx = . (2.3) 
Jn 

Assumption 1. We assume that X lies in the plane generated by ei and 62. 
Here and below we denote by ei , 62 , 63 the vectors of the canonical base of 

Denote by C the space of the v's such that (|2.3p holds. 

Remark 2.1. In principle, the space C should be an infinite dimensional func- 
tion space, so in order to discuss dynamics one should introduce a suitable norm 
in it, prove an existence and uniqueness theorem for the solutions of the Cauchy 
problem and, in order to use energy conservation (or dissipation) to prove dy- 
namical properties, one should also prove that dynamics is well posed in the 
energy space, which is in general unknown 113 j. 

In the present paper we do not want to enter such a kind of mathemati- 
cal problems, so we assume that the effective number of degrees of freedom is 
arbitrary but finite, i.e. C is finite dimensional. 

For example this can be obtained by cutting at an arbitrary order the mul- 
tipole expansion of the configuration v. 

Then we would like to factor out rotations in a way similar to translations, 
however this requires a careful discussion. The point is that it is clear what it 
means to rotate a body, but it is not clear how to say that a deformation does 
not rotate the body: as we will see, this is not a well defined concept. 

The first point is that there exist two different ways of rotating an elastic 
body. These are described by the following two actions of 5*0(3) on C 

Ai : 50(3) X C ^ C 

(r,w) ^Tv (2.4) 

and 

A2 ■■ 50(3) X C ^ C 

{R,v) RvoR-^ . (2.5) 
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The action Ai consists in rotating the elastic body rigidly in space. Clearly this 
kind of rotation is not what one would call an elastic deformation. 

The action of A2 transforms a configuration in which the body is deformed in 
a given way, into a configuration in which the deformation is the same, but in a 
different body direction. We mean that if, for example, the initial configuration 
is an ellipsoid with some principal axes, then the second one is an ellipsoid with 
the same shape, but with axes which have been rotated inside the body. This 
is a true elastic deformation that involves dissipation. 

The couple of group actions just defined allows one to introduce a structure 
of principal fibre bundle in C, the base manifold being the quotient A4 :— 
C/[S0{3) X 5*0(3)]. The elements of such a quotient manifolq^ are what one 
could call "deformations which do not rotate the body". 

As usual in this geometric context, it is useful to introduce coordinates in 
which a point of C is represented by an element of SO{3) x 50(3) and an element 
of the base manifold. However a concrete representation of the elements of the 
quotient manifold can be obtained only locally, by introducing a local section 
of the bundle, namely by choosing a submanifold of C transversal to the group 
orbit. This will be done in section U7I\ by introducing the section S. Obviously 
there are infinitely many possible choices. 

Remark 2.2. In the literature only the action (|2.4p is usually mentioned, and 
sometimes \15].\lGl the section is chosen by saying that (in the case of a homo- 
geneous body) the elements of the section satisfy 

[ cni\{v{x))d^x=0 . (2.6) 

Of course this is a possible choice. All we will do is independent of the choice 
of the section. 

Since the action A2 is not fre€0, we will consider also the action ^3, combi- 
nation of the actions Ai and A2^ defined by 

A3iR)v := AiiR)A2iR^^)v = V o R (2.7) 

and study the couple of actions Ai and ^3. 

The advantage in introducing ^3 is that such an action is free, at least in a 
neighborhood of the identical deformation 

wo(x)=x. (2.8) 

2.2 Tensor of inertia and adapted coordinates 

We recall that the matrix of inertia of the body is a symmetric matrix / = 
{lij}^ j^i whose elements are 

hj = Itj{v) • / w(x) A (e^ A v{x))p{x)d^x . (2.9) 

Jn 

^To be precise, since the group action is not free (as we will discuss below), Ai is not a 
manifold, but this is not relevant for the present discussion. 

fact which gives rise to some problems when trying to introduce coordinates 
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Tie eigenvalues of such a matrix are called the principal moments of inertia and 
will be denoted by Ii, 12,13- The eigenvectors of I are called principal axes of 
inertia and will be denoted by Ui , U2 , U3 . 

We introduce now an adapted set of coordinates in a neighbourhood of the 
"identical" deformation vq (excluding however such a configuration). To this 
end we need to introduce a few objects: 

(1) Define 

C^ = {ve C\h ^ hJi ^ I3J2 ^ h} (2.10) 
and its complement 

C= = {v e C\l, = Ij for some i ^ j) . (2.11) 

This is useful since the principal axes Ui, U2, U3 are uniquely determined 
in . 

(2) Define 

V := {v e C\I{v) is diagonal} . (2.12) 

We also define T>^ := T> DC^. Observe that T> is a codimension 3 sub- 
manifold of C, invariant under the action A3 (we will show in the proof of 
Lemma 12.81 that the action A3 leaves invariant the matrix of inertia) and 
observe that vq ^V. Moreover, as we will prove in Lemma [^751 onTX^C^ 
the action Ai is independent of the action ^3, and is transversal to V. 

(3) Consider the group orbit yl3(50(3))uo C I?, and let 5 C I? be a codimen- 
sion 3 (in I?) manifold transversal to such a group orbit. Actually we are 
interested in the restriction of such a section to a small neighbourhood of 
■yo- We still denote by S such a local section. The existence of such an 
S is assured by the fact that the action A3 is free and therefore defines a 
foliation of T). 

(4) Finally define F to be the tube constituted by the orbits of A10A2 starting 
in 5 n , namely 

J- Ai{SO{3))A2{SO{3)){S n C^) . (2.13) 

We remark that T = Ai{SO{3))A3{SO{3)){S D C^). 

Remark 2.3. The eigenvalues Ij, as functions of the matrix elements {hj} 
(and therefore of the configuration v), are smooth functions on C^; however, 
the first derivatives of the Ij 's have a singularity at C=, therefore the Ij 's can 
be used as Lagrangian coordinates only on . 

Remark 2.4. When restricting to the submanifold T> , the eigenvalues Ij co- 
incide with the matrix elements on the main diagonal, and therefore they are 
obviously smooth functions of the configuration. 

Assumption 2. We assume that the functions Ij : T) ^ M, j = 1,2,3 are 
independent in a neighbourhood of vq . 

Remark 2.5. This is true if, for instance, for any j there exists a deformation 
which changes Ij, but does not change Ii, i ^ j . For example, this could be false 
in the incompressible case. 
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In the rest of the section we will prove the following Theorem 
Theorem 2.6. There exist functions (zi, Z2, ■ ■ ■), with 

zy.S^R (j = l,2,...) 

such that: 

(i) {Ii, I2, 13, zi, Z2, ■ ■ ■) is a set of smooth coordinates on S. 
(ii) the map 

50(3) X 50(3) X (5nC^) 9 (r, i?, /i, /2, h,zi,z2, . . .) ^ Ai{T)A2{R)w e F 

(2.14) 

is a 24-fold covering of J- ; here we denoted w — {Ii, I2, 13, zi, Z2, . . .). 

Remark 2.7. The number 24 arises as the order of the chiral octahedral group. 
More precisely it corresponds to the number of ways in which an oriented triple 
of orthonormal vectors can he rotated in such a way that the vectors lie on a 
triple of unoriented fixed orthogonal axes. 

2.2.1 Proof of Theorem [2761 

To begin with, we prove the existence of (21, Z2, . . .) satisfying (i). Observe that 
iS is a smooth submanifold of C. Then, since Ii{w), /2(w), hiw) are independent 
functions, it is possible to complete the triple (Ii, 12,13) to a local system of 
coordinates (Ii, I2, 13, zi, Z2, . . .) near vq. 

In the rest of the section, we will prove (ii). As a first step, we show that the 
two actions of 50(3) on are independent, which is implied by the following 
Lemma. 

Lemma 2.8. For any fixed v G C^, we consider two subspaces ofTyC^, tangent 
to the group orbits Ai{SO{i))v and A3{SO{3))v , namely 

Ti := nAiiS0{3))v T3 := T«^3(50(3)){; . 

Then, 71 H Ti = {0}. Moreover, if v €z D T>, then 7i is transversal to T> . 

Proof. In order to prove the thesis, we start by showing that the action Ai 
rotates the matrix of inertia, while the action ^3 leaves it invariant. We have 

I,j{Tv) = e, • / ru(x) A (e^ Ar{;(x))p(x)d^x 
Jn 

= r-^ei • / w(x) A (r^^ej Au(x))p(x)d^x , (2.15) 
Jn 

which shows that the I{Ai{T)v) is the matrix of I{v), just referred to a rotated 
basis. On the other hand, we have 

Uj{v oR)^e,- / v{Rx) A (e^ A u(i?x))p(x)(i^x . (2.16) 
Jn 

Setting y = i?x, we have 

h,{v o i?) = e, • / v{y) A (e, A f;(y))p(y)d3y , (2.17) 
Jn 
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which means that the action of A3 leaves the matrix of inertia invariant. This 
imphes 

dl{v)v3 = Vw3 e Ts , (2.18) 

while 

dl{v)vi ^ Vvi e Ti \ {0} , (2.19) 

from which the independence follows. 

To get the transversality when v ^ O V, we remark that Ai rotates the 
principal axes of inertia, then, since the three eigenvalues are distinct, it destroys 
the diagonal structure of /. □ 

Remark 2.9. As an obvious corollary of Lemma \2.fA we also have that Ai and 
A2 are independent at any point v S C^, in the sense that 7i is transversal to 
the tangent space T2 '■= T,uA2{S0{'i))v . 

Moreover, we observe that in the three eigenvalues of the matrix of inertia 
are distinct, so the eigenvectors Ui, U2, U3 are well determined. Furthermore 
the dependence of the eigenvalues and eigenvectors on the configuration v G 
is smooth. 

Now, we want to represent any configuration v £ F m. the form 

v = Ai {V)A2 {R)w , weSC^C-L. (2.20) 

Let us first represent any 1; G in the form 

V = Ai{f)w , w G . (2.21) 

Proposition 2.10. The map 

n : SO{i) xV^^ 

it,w)^ Ai{f)w 

is a 24-fold covering map. 

The proof will make use of the following Theorem, which is an immediate 
corollary of 0, Proposition 1.40, p. 72. 

Theorem 2.11. If G is a finite group, acting freely on a Hausdorff space X , 
then the quotient map X — > X/G is a covering map. 

Proof of Proposition \2.1(A We observe that each w G has many distinct rep- 
resentations of the form (j2.2ip : since the three principal moments of inertia 
are distinct from one another, the directions of the principal axes of inertia are 
well-determined, but the same is not true for what concerns their orientation; 
moreover, any of the principal axes may be labeled Ui as well as U2 or U3. In 
order to make this rigorous, consider the equation 

Ai{fi)wi^Ai{f2)w2, (2.22) 

with wi,W2 € ■ This implies 

^2 =^i(f^^fi)t&i . (2.23) 
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Therefore, since wi,W2 G 2?^, the rotation f^'^fi must transform the set 
{ui,U2,U3} to a set of miit vectors having the same directions. It is easy 
to see that the set of rotations satisfying this property is the subgroup of 5*0(3) 
generated by the three rotations 
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Such a subgroup, which we wih denote by O, is isomorphic to the group of the 
orientation preserving symmetries of the cube, which is a group of order 24, 
known as the chiral octahedral group. This argument shows that the possible 
representations of the form (|2.2ip are at most 24. On the other hand, for any 
w E T)^ and F G 5*0(3), we have that the expression 



A(fro)[^i(roi)7i] 



(2.24) 



yields 24 different representations of the same configuration, as Vq varies within 
the group O. Therefore, each configuration v e has exactly 24 distinct 
representations of the form (j2.2ip and a natural identification arises between 
and (50(3) x V^jO, where the action of O on 50(3) x is defined by 



[ro,(f,zi)]^(fro,A(roi)w) . 



(2.25) 



Now, applying Theorem 12.111 with X 
thesis. 



50(3) X and G = O, we get the 

□ 



The above Proposition given as a global statement applies also to a small 
tube of orbits originating in S. Precisely, define T '■— y^3(50(3))(5 nC^): then 
we have 



Corollary 2.12. 

n :50(3) X r 

(f,^i) 

is a 24-fold covering map. 



Proof. The only thing we have to prove is that is the image of 50(3) x T 
through n. However, this is obvious, since 



n(50(3) X r) - AiiS0{3)){T) = Ai{SOi3))A3{SO{3)){SnC^) = T 



(2.26) 
□ 
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End of proof of Theorem \2.6\ The last step consists in factoring out the group 
action ^3. This is easy, since the action is free. Therefore, one can decom- 
pose 

T3w=A3{R)w (weSnC^) (2.27) 
in a unique way, and moreover the map 

w I— > (R, w) 

is smooth. Therefore, a 24-fold covering of T is naturally induced by the map 

(f, R, w) ^ Aiit)A3iR)w . (2.28) 

Now, setting 

r := 
R :~ R ^ , 

we find that also 

{T,R,w)^ AiiT)A2iR)w (2.29) 
is a 24- fold covering of J^, which completes the proof of (ii) and of Theorem l2.6l 

2.3 Elastic potential energy 

Let us study the form of the elastic potential energy and of the potential en- 
ergy of self-gravitation in the coordinates just introduced. With an abuse of 
terminology, we will call the sum of these two potential energies simply "elas- 
tic potential energy" and we will denote it by T4; furthermore, we will refer to 
the corresponding forces as to the "elastic forces", leaving understood that they 
include also the forces related to self-gravitation. 

Assumption 3. The identical deformation Vo{x) = x is a minimum of the 
elastic (and self- gravitational) potential energy. 

Remark 2.13. This Assumption means that we have chosen the material space 
51 as representing an abstract realization of the equilibrium configuration of the 
compressed elastic sphere under the effect of self- gravitation (and not of the 
totally undeformed elastic body). 

In the equilibrium state, because of the rotational invariance, all three prin- 
cipal moments of inertia are equal to the same constant Iq. For simplicity, we 
use the differences between the Ij 's and /q as configuration variables instead of 
the Jj's themselves, so we define 

J,:-/, -Jo, (i = 1,2,3), (2.30) 

and we assume (without loss of generality) that 

Zj{vo) =0 V j . 

Remark 2.14. Due to the Ai- and A2-invariance, the elastic potential energy 
does not depend on T and R. 

We also assume that the minimum is nondegenerate and that the body is 
very stiff. Summarizing we make the following Assumption: 
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Assumption 4. The elastic potential energy has the form 

K( J, z) = -Vo{J, z)=- [Q(J, z) + VsiJ, z)] , (2.31) 

e £ 

where e is a small parameter, Q a nondegenerate quadratic form and V3 has a 
zero 0/ order three at the origin. 

We want to study more in detail the form of the elastic potential near the 
equilibrium, but we have to cope with the fact that our coordinates are singular 
at the equilibrium configuration wo(x) = x. 



Lemma 2.15. The elastic potential energy has the form 

A 
2 



QiJ,z) = ^{Jl^ + J2^ + J3^) + B{JiJ2 + JlJ3 + J2J3)+ (2.32) 



+ ^CjZj(Ji + J2 + J3) + 2 X! DjkZjZk , 
j=i j,k=i 

where the constants A, B,Cj, Djk are such that the quadratic part Q{J,z) is a 
positive definite quadratic form in the variables (J, z). In particular, this implies 
A> B. 



Remark 2.16. By Remark \2JA such an expression can he used to compute the 
Lagrange equations only outside C= . 

Proof. Any v ^ J- can be represented as Ai{r)A2{R)'w, for some T, R G 5*0(3) 
and w € S D C^. Moreover, due to the group action invariance, the potential 
energy associated to the configuration v must be the same as the potential 
energy associated to the configuration w. Therefore the functional form of the 
elastic potential in terms of the variables J, z can be computed working in S 
and then the obtained form holds on the whole of J^. 

The independence of the Lagrangian of the body with respect to the choice 
of the representative (F, i?, w) from the covering mentioned in Proposition 12.61 
implies that the expression of the elastic potential energy must be symmetric 
with respect to permutations of the indices i = 1,2,3, and the expression of 
Q{J, z) in (j2.32p is the most general expression of a quadratic form with such a 
property. □ 

2.4 Planar restriction 

We are going to study the dynamics of the satellite in the special case when 
the spin axis of the body is orthogonal to the plane of the orbit and coincides 
with one of the principal axes of inertia of the body. For this reason, in the rest 
of the paper we will restrict to "planar deformations" and "planar rotations". 
Precisely, we make the following assumptions. 

Assumption 5 (Planar deformations). The configuration is such that 63 is an 
eigenvector of I . We label the principal axes of inertia so that 1*3 63. 

In other words, we are assuming the matrix of inertia of the satellite to be 
of the form 

hliv) Il2{v) 
Hv) = h2{v) l22{v) , (2.33) 

hsiv) 
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so that hiv) = hsiv). 



Assumption 6 (Planar rotations). F is a rotation about the e^-axis, i.e. it has 
the form 

cos a — sin a 
r = r(a):= sin a cos a . (2.34) 
1 

As a consequence of these two assumptions, R is also a rotation about the 
ea-axis, i.e. there exists an angle (3 such that 



R = i?(/3) 



cos /3 — sin (3 
sin /3 cos (3 
1 



(2.35) 



Remark 2.17. The assumptions \^ and \^ together with Theorem \2.b\ imply 
that (a, /3, /i, I2, 13, Zi, Z2, ■ ■ .) are good Lagrangian coordinates for the body con- 
figuration. Actually, by following the proof of Theorem \ 2.b\ one can show that 
.such coordinates form a 4-fold covering of the configuration space restricted to 
planar configurations. 

The fact that dynamics remains confined for all times within the set T will 
be guaranteed by the local stability result proved in the following sections. 



2.5 Kinetic energy 

By Konig's second Theorem, the kinetic energy T can be written as the sum of 
two terms: the former is the kinetic energy of the center of mass 



J- r.l 



(R^ + i?>^) 



(2.36) 



and the latter is the kinetic energy of the satellite with respect to its center of 
mass 

Tr := Tr{a, 13, Ji, Js, z; Ji, J2, J3, z) . (2.37) 

Remark 2.18. T^ is independent of a and /3 due to the rotational invariance 
of the satellite. 



We will use the notation 



^ 5+n 

:= - ^ aik{J,z)qiqk , 



(2.38) 



i,k=l 



where q = (a, /?, Ji, J2, J3, z). Observe that the coefficients aik{J,z) are such 
that the quadratic form is positive definite on T . 



Lemma 2.19. The coefficients aik{J,z) satisfy: 
(i) 



(ii) 



aii(J, z) = h = I0 + 



lim ai2( J, z) = 

(J,z)->0 



(2.39) 
(2.40) 
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Proof. For w G S, we set 

Wde/(x) := w(x) - t;o(x) = w(x) - x . 

and 

U ^ A2{R{fi))Wdef ■ 

We remark that A2{R{l3))vo = uq. 

Now, let us evaluate the kinetic energy Tr. For v Cz J- we have 

V = r(a)(x + ii(x)) . 

Taking the derivative with respect to time, we get 

w(x) = — ^(x + u(x)) + r(a)u(x) . 



(2.41) 
(2.42) 

(2.43) 
(2.44) 



Therefore, 

Tr - 



\ [ [^(x)] Vx)d^ = i / [r(-a)t;(x)]V(x)d3x: 



T{-a) 



dVia) 
dt 



(x + w(x)) + 'u(x) 



/9(x)d^x 



{w X [x + u(x)]}^ p(x)d^x + 



+ / (a; X [x + u(x)] ,u(x))p(x)d3x+ - / [?i(x)] p(x)d^x,(2.45) 

where u> is the angular velocity of the satellite, defined by wx (•) = [r(— a)] [■^T{a)]{). 
Under our assumptions, lo — du3. 

As the vector field u(x) is independent of a, we observe that Tr is the sum 
of three integrals, the first of which gives the term in 6? , while the third one is 
a quadratic form in (/3, Ji, J2, J3, i) and the second one gives mixed terms in a 
and in the other velocities. 

Therefore, one gets 



aii{J,z)= I {U3 X [x + M(x)]}^p(x)rf^x , 
Jn 



(2.46) 



and it can easily be seen that this expression equals the moment of inertia 
related to the vertical axis, thus proving (i). 

For the second part of the Lemma, we have to study the coefhcient of the 
term d/?. Observe that such a term arises from the integral 



(a; X [x + u(x)] , 'it(x))p(x)d^x 



in (j2.45p . Here, the a factor comes from the angular velocity w, while the /3 
factor is hidden in m(x). Using 



u{^) = R[f3)wdef[R{-f3)^] 



(2.47) 
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we get 

^x(x) = R{l3){waef[R{-/3)^]} + $^^Wdef[R{-l3)^] + 

+pR{0)Vwaef [i?(-^)x] • ^^l^x . (2.48) 

Here, we notice that the first of the three addenda is independent of $, so we 
have 



ai2{J,z) 



i ^ (U3 X {x + i?(/3)«;,e/[i?(-/3)x]}, ^^w,ef [i?(-/3)x] + 
+ R{l3)Vwaef [i?(-/3)x] • ^^Mx)p(x)d3x , (2.49) 



which goes to zero when Wdef 0, i.e. when {J,z) 0; this completes the 
proof of (ii). □ 

3 The Lagrangian of the system 

In this section wc write the Lagrangian of the complete dynamical system con- 
sisting of the satellite subject to the gravitational field of a pointlike center 
(planet) having mass M. 

3.1 Gravitational potential energy 

To start with, we fix some notation: (R, ip) are the polar coordinates of the 
center of mass of the satellite in the plane of the orbit. We denote with 7 the 
angle between the principal axis Ui and the line joining the planet to the center 
of mass of the satellite. Such a line is usually referred to as the line of centres. 
Observe that 7 = a + /3 — ■(/;. Furthermore, we set x := a — f/', i.e. x is the angle 
of rigid rotation of the satellite, measured with respect to the line of centres. 

Proposition 3.1. In the quadrupole approximation the gravitational potential 
energy of the body in the field generated by the mass M is given by 

C Mm C M 

Vg{X, Ji, J2, J3, 7) := - + [-Ji + 2 J2 - J3 + 3( Ji - J2) cos2 7] (3.1) 
Proof Define p : ({n) — ^ M as 



detli(C-HC)) 



i.e. p(^) is the density of the satellite at the point ^ = C(x). Let {xi,X2,X3) be 
the Cartesian coordinates referred to the system with origin X and axes U1U2U3. 
Then we introduce the spherical coordinates (r, 1?, 0) of the generic point P in 
the satellite, defined by: 

xi = r cost? cos (/) (3.2) 
X2 = rsiwd coscf) (3.3) 
Xz = rsin</) . (3.4) 
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In this frame, the products of inertia lij vanish, i.e. 

/ p{^)r^ cos^ (l)COS-dsmd(f^ = (3.5) 
Jan) 

/ pi^r"^ cos (l> sin (j) cos 'ffcf^ = (3.6) 

Jan) 

/ p(Or^cos0sin(/)sinM^^ = 0, (3.7) 
Jan) 



'an) 

and the principal moments of inertia are given by 



h = Io + Ji^K2 + K3 (3.8) 
I2 = Io + J2 = Ki+K3 (3.9) 
^3 = Io + J3 = Ki+K2, (3.10) 



where 



Ki = I- [ p{^ycos'^(t>cos'^M^^ (3.11) 
2 Jan) 

K2 = I [ p(0r^cosVsin2-i?d3^ (3.12) 
2 Jan) 

-- I [ Pi^Vsin'cj^d'i. (3.13) 



K: 



The gravitational potential energy Vg is: 

V, = -[ ^dH = -f ^^M^=,^^, (3.14) 

Jan) I? I J(;{n) y/R^ + r'^ -2 Rr COST] 

where rj is the angle between the line of centres and XP. Notice that the relation 

cos r/ = cos cos (1? + 7) (3.15) 
holds. Let us recall now how the multipole expansion arises. We have 



lel y/R^+r^-2Rrcosri R + (^^^^ _ 2 cosr, 
In terms of the Legendre polynomials Pniz), one has 

^a;"P„(^). (3.17) 



In particular, we recall that 

Poiz) = 1 
Pi{z) = z 
3z^-l 



P2{Z) = 
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Taking the quadrupole approximation means to cut the sum at n = 2. We get 

1 1 ^r- , , 1 



n>0 

SO the potential energy becomes 



-V^'^^'^-^' , (3.18) 



C(n) ^ 



Here, the first term equals — the second one vanishes because X is the 

center of mass of the satelhte; the third term, namely 

[ GMp(0(3cos^r?-l)r^ 3, 
gives what we call the "tidal" potential energy. A brief manipulation shows that 
-S^ Jan) 2 



i?3 



/ p(^)r^[3cos^ <?!)(cos 1? cos 7 — sin'!?sin7)^ — l\d^^ = 

Jc(n) 

r CM f 

-{K, cos' J + K2 sin^ + piOr'd'^ = 

Jan) 



■ „o [Ki cos 7 + A2 sm 7) 



+ ^ / p(0^'[sin^<^ + cos2(/)(sin2^? + cos2t?)]c/3^ = 

Jan) 

= ^(^1 I + K2 sm' 7) + -^(^1 + K2 + K3) = 

= 2^2 -J3 + 3( Ji - J2) cos^ 7] • (3.20) 

□ 

Remark 3.2. This result is probably well known, but we did not find a reference 
for it. 

4 Conservative dynamics 

After having described our model of satellite, we are ready to study the dynamics 
of our system. At first, we neglect the dissipative effects and study the dynamics 
of the corresponding conservative system. 
The Lagrangian 

C = T-V (4.1) 

of the system is the difference between the kinetic energy T and the total po- 
tential energy. Therefore, collecting the above results one has: 

_ m 
^ ~ 2 



+ R'lP'') +Tr{x + i^, j, /3, i; J, z) + (4.2) 



GMm GM 



+ — ^ + ^rT [-^1 - 2 J2 + J3 + 3( J2 - Ji) cos2(x + /?)] - V,{J, z; e) . 
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Now, we observe that the Lagrangian does not depend on the cycHc coordinate 
■0, so the total angular momentum 



dC 



5+n 



p:= -J =mR'^^ + iJ3 + Io)ix + i^) + 2y2aik{J,z)qk , (4.3) 



k=2 



is a constant of motion. We can invert relation (|4.3p . to get the expression of V' 
as a function of the other variables: 



mR'^ + Ia + J3 

Then, we can drop one degree of freedom and study the reduced Lagrangian 



(4.4) 



(4.5) 



where ^ must be thought of as a function of the other Lagrangian coordinates 
and velocities. After some calculations, we get 



£* ^T2 + Ti-V 



where 



T2 - 

V = 



-R^ + Tr{x,jJ,i;J,z) 



(^3 + ^)x + 2 J2kt2 aik{J, z)qk 



{J3 + Io)x + 2 o-ik{J, z)qk 



2(mi?2+/o + J3) 



(4.6) 



-(4.7) 



mi?2 + /(, + J3 



P 



2(mi?2 + I0 + J3 
GMm GM 



(4.8) 
(4.9) 



R i?3 
Such a system has the conserved quantity 

5+n 



[Ji - 2J2 + J3 + 3( J2 - Ji) cos2 7] + Ve{J, z; e) 



dC* 



k=l 



(4.10) 



where 



y {R, Xi Ji, J2, J3, zi,z2, ■ . .) 



and the strict minima of V are Lyapunov-stable equilibria of the system. 
Let i?o be a nondegenerate minimum of the function 



VaoiR) 



GMm 



P 



R 2{mR^ + h) ■ 
Then we have the following 

Lemma 4.1. For any small enough e there exist R, J,z, s.t. 



(4.11) 
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(1) the manifold 

M := {(i?,X,/3, Ji, J2,J3,^)|x + /? = 0} , 

is composed by critical points ofV. 

(2) A4 is a minimum ofV which is nondegenerate in the transversal direction. 

(3) One has (Ji, J2, J3, z) = 0(e) and\R- Rq\ ^ 0{e). 

(4) Finally Ji < J2 < J3. 

Remark 4.2. Point (iv) guarantees that M C C^. If e is sufficiently small, 
then we have Ai d J-. 

Remark 4.3. Ai is the manifold corresponding to 1:1 spin orbit resonance. 

Proof. We look for a minimum of V in the domain Ji < J2 and \ J\ < Ce for 
some fixed C. 

First remark that, as a function of 7 = % + /3, V has a minimum at 7=0 
(strict if Ji < J2). Consider now ^^|^_g; as a function of R it has a nondegenrate 
minimum at some point R = R{J, z) fulfilling 

|i?(J,z) - i?o| < Ce . 

Consider now the restriction V = V{J, z) oi V to the mainfold j — 0, R — 
R{J, z); since 

17( J, z) = i [Q( J, z) + VsiJ, z)] + 0(1) , (4.12) 

such a function has a nondegenerate minimum close to zero. 

Then (1), (2) and (3) follow provided one shows that Ji < J2. We are now 
going to prove (4) which in particular implies the thesis. 

To this end, observe that at the critical point one has 







dV 2GM A - B 



-Ji + - (J2 + J3) + - 5Z Q% + 0{e) (4.13) 



dJi R^ e e e 



r)\/ OAf A - R - - 1^ 
= TTT = + 7'^^ + -i'h + J3) + -Y. + 0{e) ■ (4.14) 



9J2 R^ e e " e 



dV _ p^ GM 

^ "2(mi?2 + /o + J3)2 " -rT + 

A - R - - 1 " 
+ - J3 + -{Ji + J2) + - ^ CjZj + 0(e) . (4.15) 



Subtracting (|4.14p from (|4.13p . we obtain 
3GAf A-B 



i?3 



(Ji - J2) + 0(e) = . (4.16) 
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The positive definiteness of the quadratic form Q imphes A — B > Q. Therefore, 
if e is sufhciently small, we have Ji < J2. Subtracting (|4.15ll from (|4.14p . we 
get 

^2 l\ ^ -J ^2 + ^(-^2 - J3) + 0(£) = . (4.17) 
2(mR^ + I0 + J3) £ 

Hence, for sufficiently small e, we have J2 < J3. □ 

Corollary 4.4. The critical manifold Ai is an orbitally stable equilibrium for 
the Lagrangian system of equations 



d dC* dC* 



(fc = l,2,...) (4.18) 



dt dyk dyk 

5 Dissipative dynamics 

In order to take into account the dissipative effects related to the internal friction 
of the satellite, we introduce a Rayleigh's dissipation function F{y;y) (namely 
we assume that the dissipative equations of motion have the form (|5.ip below) . 
Of course, we do not expect dissipation to act directly on the orbital variables 
R and Xi so we assume F to be independent of R, x, R and Due to the 
y^2-iuvariance of the satellite, F is also independent of the coordinate (3 (but it 
depends on $); however, for the proof of our result, it is not necessary that F 
be /3-independent. 

Assumption 7. As a function of the velocities {f3,J,z), the Rayleigh's dissi- 
pation function F(/3, J, i; J, z) has a nondegenerate minimum at 0. 

We denote by y'^ := (/3, J, z) the variables fixing the configuration of the 
elastic satellite. So in particular we have F — F{y'^,y'^) with dF/d/3 = 0. 
Now we can state our main result: 

Theorem 5.1. If e is sufficiently small, the manifold Ai, defined in the previ- 
ous section, is an asymptotically stable equilibrium for the dynamical system of 
equations 

dt dyk dyk dyk 

In order to prove this Theorem, we will use a common tool in the study 
of dynamical systems, known as LaSalle's invariance principle (see |11) . §2.6), 
which allows one to prove results of asymptotical stability in presence of a 
Lyapunov function E satisfying a nonstrict inequality of the type < 0. To 
state this principle, we first recall the definition of an invariant set. 

Definition 5.2. A subset M of the phase space is called (positively) invariant 
if all solutions starting in M remain in M for all future times. 

We now state a version of LaSalle's Theorem. 

Theorem 5.3 (LaSalle's invariance principle). Suppose that E is a real-valued 
smooth function defined on the phase space, satisfying E{y, y) < for all {y, y). 
Let M be the largest invariant set contained in |(?/, y)\E(y, y) ^ o|- Then every 
solution that remains bounded for i > approaches M as t ^ 0. 
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Proof of Theorem \5.1l Let us evaluate the time derivative of the energy: 

1^ - —(\^- —\ - —r* - 
Itt ~ di [f^^^'^d^ J ' di 

^ dc* ddc* dc* ^ dc* 

= -Ey'^l^--(^Q^ + T.y^p-] ' (5-2) 

where is the quadratic part of F, and F3 is the part of order 3 in y"^. Let 
J\fV := {{y,v)\F{y;y) = 0} = y)|/3 = j = i = o| be the subset of phase 
space where there is no energy dissipation. Then, due to LaSaUe's invariance 
principle, any solution such that (i/(0), 2/(0)) belongs to a sufficiently small neigh- 
borhood of (A^ , 0) (notice that such a solution will stay bounded for alH > due 
to the Lyapunov stability of M which has been proved in the previous section) 
will get arbitrarily close to the largest invariant subset of AfT), for t +00. 
Therefore, the only thing we have to check is that the set AfT> contains no orbit, 
apart from the points of the manifold A4. To check this, observe that, if such an 
orbit existed, it would satisfy equations (|5.ip . In particular, the orbit satisfies 

d dC* dC* 

and 



d dC* dC* OF 



dt dp 9/3 dl3 ■ 
When restricting to A/"!?, these two equations become, respectively, 

{lo + h?X , 2mi?(/o + hfxR , 



ii?2 + /o + J3 (mi?2 + /o + J3: 



2 



(5.4) 



(5.5) 



and 



^,r^^^■■ 2pmi?(/o + J^)R dV 
(mi?2 + /q + Jg)^ c*7 



ai2(J, z)(/o + J3)x , 2mRai2{J, z){Io + J3)xR 



mR^ +IQ + J3 {mR^ + /q + -h) 



2 



(5.6) 



, , .. 2pmRai2(J, z)R dV , , 

(mi?^ + /o + J3)^ 07 



where we took into account that, by Assumption [71 dF/d/3 vanishes on A4. 

ai2(J,z) 
I0 + J3 



Multipliying (j5.5p by "Pi'^}^'' and subtracting (j5.6p we get 



I0 + J3J 

which, by (|2.40p . implies fi^ = 0, and therefore x + /3 = 0. Then, we have 

X = — /? = 0, since /? = on AfT). Now , substituting x = x = into equations 
and (|5.6p . we find R — 0. Finally, we observe that now we have 

^ ^ R = j = z = , (5.7) 
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which is true on the equihbrium manifold A4 only. 

We have thus proved that the only orbits contained in AfV are the points of 
the 1:1 spin-orbit resonance manifold M., which implies the asymptotic stability 
ofM. □ 
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